Multicanonical Study of the 3D Ising Spin Glass l 
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We simulated the Edwards-Anderson Ising spin glass model in three dimensions via the 
recently proposed multicanonical ensemble. Physical quantities such as energy density, spe- 
cific heat and entropy are evaluated at all temperatures. We studied their finite size scaling, 
as well as the zero temperature limit to explore the ground state properties. 
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One of the great challenges in numerical simulations is the study of spin glasses. These 
systems display a rather peculiar structure due to the effects of disorder and frustration, 
for reviews see [JTj, 0]. Recently two of the authors explored the multicanonical ensemble 
H for simulations of spin glasses. Results for the 2D Edwards- Anderson Ising spin glass 
were encouraging ||. Similar concepts were also tested in f|. In this letter we present 
a preliminary multicanonical simulation of the physically more interesting 3D Edwards- 
Anderson Ising spin glass. 

The Hamiltonian of the model is given by 

H = - J2 J t jSiSj, (1) 

where the sum goes over nearest neighbours and the exchange interactions Jij = ±1 between 
the spins Sj = ±1 are randomly distributed over the lattice with the constraint J2 Jij — 
for each realization. Despite its simplicity the model is supposed to be sufficiently realistic. 
Recent simulations of the 3D model in a magnetic field favour the mean field picture 
rather than the alternative droplet model |7|, |8| . However, it can be argued that equilibrium 
at sufficiently low temperatures has not been reached ||. For previous simulations of Ising 
spin glasses in 3D see JTT], [II], [12], [IB], |TJ|. 



The multicanonical ensemble || can be defined by weight factors 

P M (E) = exp[~(3(E)E + a(E)]. (2) 



a(E) and (3(E) are to be determined such that for the chosen energy range E min < E < E 
the resulting multicanonical probability density is approximately flat: 



max 



Pmu(E) = c mu n{E)P M {E) « const (3) 

where n(E) is the spectral density. In the present study we take £ max = {(3(E) = for E > 
E max ) and E min = E° the ground state energy of the considered spin glass realization. The 
purpose of the function a(E) is to give (3(E) -1 the interpretation of an effective temperature. 
This leads to the recursion relation 

a(E - 4) = a(E) + {(3{E - 4) - (3{E)} E, a(E max ) = 0. (4) 

The multicanonical function (3(E) is obtained via recursive multicanonical calculations. One 
performs simulations with j3 n (E), n = 0, 1,2, ... , which yield probability densities P n (E) 
with medians E^ edian . We start off with n = and /3°(E) = 0. The recursion from n to 

n + 1 reads 

f p»(E) for E > E n median] 



(3 n+ \E) = I 



(3 n (E) +0.25x \n[P n (E + A)/P n (E)} 

for F, n > F, > F, n ■ ^ ' 

101 ^median > ^ dL -^ m in) 



I /T +1 (£ mm )ior E<E^ m . 
The recursion is stopped for m with E™^ = E° being groundstate. 



By weighting with exp[— J3E + /3(E) E — oc(E)] canonical expectation values 0(/3) = 
Z0)~ l He 0{E)n{E) exp(-/3£), where Z((3) = J2 E n(E) exp(-J3E) is the partition func- 
tion, can be reconstructed for all J3 (the canonical inverse temperature). The normalization 
constant c mu in equation (3) follows from Z(0) = Y^,E n {E) = 2 N , where N is the total 
number of spin variables. This gives the spectral density and allows to calculate the free 
energy as well as the entropy. 

We simulated the Edwards-Anderson Ising spin glass model on 3D cubic lattices with 
linear size L = 4, 6, 8 and 12. We performed simulations of 2 x 10 6 iterations on 4 3 lattice 
with 32 different realizations of random variables Jy. We carried out up to 4 x 10 6 iterations 
on 4 3 , 6 3 and 12 3 lattices with 16, 8 and 4 realizations, respectively. Thermal averages 
were evaluated after the first additional 2 x 10 5 iterations being discarded, although with a 
disordered starting configuration the multicanonical ensemble is immediately in equilibrium. 

Our results for physical quantities are summarized in Table 1. The final mean values 
and their error bars are obtained by combining the results from the different realizations. 
Different realizations are statistically independent and enter with equal weights. The final 
error bars are enlarged by a Student multiplicative factor, such that the probability content 
of two standard deviations is Gaussian (95.5%). 





L=--A 


L = 6 


L = 8 


L = 12 


Pmax 


1.04 ± 0.03 


1.52 ± 0.12 


2.04 ± 0.18 


2.30 ±0.38 


rl 


685 ± 95 


30166 ± 16383 


117038 ± 40630 


1580230 ± 882686 


e u 


-1.7403 ± 0.0114 


-1.7741 ± 0.0074 


-1.7822 ± 0.0081 


-1.7843 ± 0.0030 


s° 


0.0724 ± 0.0047 


0.0489 ± 0.0049 


0.0459 ± 0.0030 


0.0491 ±0.0023 





0.65 ± 0.05 


0.71 ± 0.05 


0.56 ± 0.10 


0.68 ± 0.10 



Table 1: 



The energy density e(/3), the specific heat c((3), and the entropy per spin s((3) follow 
in a straightforward manner by constructing the canonical ensemble. Fig. 1 depicts the 
energy density and the entropy per spin versus inverse temperature for L = 8. The indicated 
error bars are with respect to the eight different realizations. For the temperature range 
0.8 < T < 1.6 we collect in Fig. 2 our L = 4, 6, 8 and 12 results for the specific heat. This 
is of interest as a spin glass phase transition is claimed to take place at T c « 1.2 [110, ITT 



However, this picture was recently questioned by Bhanot and Lacki |T5j on the basis of exact 
partition function calculations on fairly small systems (up to 5 x 5 x 9). The critical exponent 
estimate v = 1.3 [|1(], [TTJ implies a = —1.9 via the hyperscaling relation Dv = 2 — a. The 
finite size scaling (FSS) of the specific heat T = T c is c ~ L*/". Clearly, Fig. 2 is inconsistent 
with this L-dependence and exhibits instead almost negligible finite size effects as would be 
typical for non-critical behaviour. 

Let us now concentrate on the zero temperature (ft — > oo) limit. The groundstate energy 
density is e° = E°/N, and we obtain its entropy per spin as s° = S°/N = ln[n(E°)]/N. FSS 
fits of the form /£ = /£, ± c/N are used to estimate the infinite volume groundstate energy 
density and entropy per spin. Fig. 3 displays these fits. Our groundstate energy density 



result e 



o 



16 and e 



= —1.7863 ± 0.0028 is consistent with the previous estimates e 
11 - —1.75 |14||. Our value for the groundstate entropy per spin s° 



-1.76 ± 0.02 
0.046 ±0.002 



Figure 1: Energy density and entropy per spin versus /3 (from the L = 8 lattices). 



Figure 2: Specific heat versus temperature for L = 4(x), 6(0), 8(D) and 12(*) lattices. 



is consistent with s° = 0.04 ± 0.01 given by Morgenstern and Binder |TB[ and lower than 



Kirkpatrick's value s° = 0.062 which was estimated from a single sample of 20 3 spins \\ I 
Our groundstate entropy value translates, even for moderately sized systems, into large 
numbers of distinct groundstates. For instance s° = 0.046 implies approximately 3.3 x 
10 34 groundstates for a 12 3 lattice. We also tried the entropy fit Ns° L = ln(iV) + c, which 
corresponds to a power law ansatz for the groundstate degeneracy. Unacceptable x 2 values 
rule out this fit even if the lattice range is restricted to L = 4 — 8. 

To visualize the low temperature behaviour, we show in Fig. 4 the spin glass order 
parameter distribution of one of our L = 8 realizations. Its calculation follows the lines of 
M. One clearly sees five configuration space valleys which are separated by high tunneling 
barriers. The multicanonical simulation overcomes these energy barriers by connecting back 
to the disordered high temperature states. The realization from which Fig. 4 is depicted 
showed nine tunneling events. Here tunneling means that starting from the high temperature 
disordered phase the algorithm finds a true ground state and then travels all the way back 
to the disordered phase. The average number of sweeps needed for this process is our 
ergodicity time r£. Table 1 shows how the ergodicity times increase with the lattice size for 
our simulations. The data corresponds in CPU time to a slowing down ~ \/ 3 ' 4< - 2 - ) which is 
similar to the one involved in our previous 2D simulations. 

Further, we have included in Table 1 our zero temperature (groundstates) results x° q 
for the spin glass susceptibility density Xq =< I 2 > /^ (note that < q >= due to the 
symmetry of the Hamiltonian) . The numbers are modestly indicative for a non-trivial Po(q) 
distribution (lack of self-averaging). We also looked at the temperature dependence of the 
spin glass susceptibility and the corresponding Binder parameter. In the neighbourhood 
of the suspected T c our results are consistent with [10,11], but our present accuracy does 
not allow a detailed finite size scaling investigation. Considerable future improvements are 
feasible. 



Figure 3: FSS estimates of the infinite volume groundstate energy density and entropy per 
spin. 



Figure 4: Spin glass order parameter distribution (L = 8 lattice). 

Let us state the conclusions. Our results on specific heat, shown in Fig. 2, raise some 
additional doubts about the standard picture of the spin glass phase transition. But on 
a qualitative level we find a bifurcation of the spin glass order parameter, as depicted in 
Fig. 4. A future, more quantitative, investigation of this quantity is desirable. Fig. 4 and 
our preliminary results for the spin glass susceptibility favour the mean field picture of spin 
glasses. To reach seminal results one has to investigate larger lattices and, first of all, to 
clarify the nature of the transition. 
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